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Abstract 

The presence of a light b squark (mi ~ 4 GeV) and gluino (nig ~ 15 GeV) might ex- 
plain the observed excess in 6-quark production at the Tevatron. Though provoca- 
tive, this model is not excluded by present data. The light supersymmetric parti- 
cles can induce large flavor-changing effects in radiative decays of B mesons. We 
analyse the decays B — > X s ^y and B — > X sg in this scenario and derive restrictive 
bounds on the flavor-changing quark-squark-gluino couplings. 
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1 Motivation 



The measured 5-quark production cross section at hadron colliders exceeds next-to- 
leading order (NLO) QCD predictions by more than a factor of two. While it is con- 
ceivable that this discrepancy is due to higher-order corrections, the disagreement is 
surprising since NLO calculations have been reliable for other processes in this energy 
range. Berger et al. have analysed 6-quark production in the context of the Minimal 
Supersymmetric Standard Model (MSSM) and find that the excess in the cross section 
could be attributed to gluino pair-production followed by gluino decay into pairs of b 
quarks and b squarks, if both the gluino and the b squark are sufficiently light In 
order to reproduce the transverse-momentum distribution of the b quarks, the masses of 
the gluino and light 6-squark mass eigenstate should be in the range m- g = 12-16 GeV and 
mi = 2-5.5 GeV. The masses of all other supersymmetric (SUSY) particles are assumed 
to be large, of order several hundred GeV, so as to have evaded detection at LEP2. 
Interestingly, a renormalization-group analysis in the framework of the unconstrained 
MSSM shows that a light b squark is most natural if it is accompanied by a light gluino 
with mass of order 10 GeV |f2j . 

Berger et al. have further observed that a light b squark could have escaped direct 
detection. For example, the additional contribution to the e + e~ — > hadrons cross section 
at large energy would only be about 2% and hence difficult to disentangle. The pair- 
production of light scalars would alter the angular distribution of hadronic jets in e + e~ 
collisions, but the present data are not sufficiently precise to rule out the existence of this 
effect []l||. On the other hand, there are important Z-pole constraints on the parameters 
of this model. Most importantly, production of the light b squark at the Z pole has to 
be suppressed, which implies a stringent constraint on the mixing angle 9 relating the 
sbottom mass and weak eigenstates ||. More recently, several authors have studied loop 
effects of the light SUSY particles on electroweak precision measurements [§], [5|, ||, finding 
potentially large contributions to the quantity Rb- However, a conflict with existing data 
can be avoided by having some of the superpartner masses near current experimental 
bounds, or by allowing for a new CP- violating phase in the SUSY sector ||. 

The null result of a CLEO search for the semileptonic decays B — > D^lir and 
B — > D(*Hx° of sbottom hadrons implies that the branching ratios for the decays b — > cl 
induced by i?-parity violating couplings, or b — > clx° with an ultra-light neutralino x°-> 
must be highly suppressed 0. However, a light b squark would be allowed to decay 
promptly via hadronic i?-parity violating couplings in the modes b — > cq or b —>■ uq 
(with q = u,s). Alternatively, it could be long-lived, forming 6-hadrons. An interesting 
consequence of hadronic .R-parity violating decays would be the abundant production 
of light baryons. This could significantly alter the thrust-axis angular distribution for 
continuum events containing baryons at the B factories. 

A striking manifestation of the light 6-squark scenario would be the production of like- 
sign charged B mesons at hadron colliders, because the Majorana nature of the gluino 
allows for the production of bbb*b* and bbbb final states J1J. Another potential signature 
at hadron colliders is an enhanced yield of tibb events ||. It has also been pointed out 
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that sbottom pairs would be copiously produced in T(nS) — > bb* and Xbj — > bb* decays 
H [l(| . Precise measurements of bottomonium decays could lead to new constraints on 
the squark and gluino masses. 

The presence of light SUSY particles alters the running of a s , and it is often argued 
that this would exclude the existence of light gluinos. This argument is incorrect. First, 
a gluino with mass m g ~ 15 GeV would have a relatively small effect on the evolution 
of a s . Taking, for instance, a s (mb) = 0.21 (a value in agreement with all low-energy 
determinations of the QCD coupling) and including the contribution of the gluino octet 
to the P function above the scale fig = rrig yields a s {mz) = 0.126, which is about 
three standard deviations higher than the canonical value a s (mz) = 0.118 ± 0.003. 
However, considering that at leading order only virtual gluino pairs contribute to the (5 
function, a more realistic treatment would include the gluino contribution above a scale 
fi g = Im-g ~ 30 GeV, in which case a s {rn z ) = 0.121, in good agreement with the standard 
value. Secondly, it is important to realize that even a value of a s (mz) significantly above 
0.118 would not rule out the model, the reason being that the characteristic scale fi 
inherent in all determinations of a s {fi) is typically much smaller than the total energy. 
This is true, in particular, for the determinations based on event-shape variables. In 
practice, the measurements fix a s (ft) somewhere between a fraction of the Z mass down 
to several GeV, where the gluino contribution to the (3 function is negligible. Using these 
determinations to quote values of a s {mz) (as is routinely done) assumes implicitly that 
the coupling runs as predicted in the SM. Finally, a careful analysis of the running of a s 
in the presence of light SUSY particles would have to include, for each observable, the 
modifications in the theoretical formulae due to virtual and real emissions of the new 
particles. These corrections could be significant, and could partially compensate effects 
arising from the modification of the (3 function. 

If we are to take the possibility of a light b squark and light gluinos seriously, then the 
theoretical study of their impact must be extended to the phenomenology of weak decays 
of the b quark. New sources of flavor violation arise from s-b-g and d-b-g couplings. The 
overall scale of SUSY flavor- changing interactions originating from gluino exchange is set 
by the factor g 2 s /m^ which is much larger than the corresponding factor Gp ~ 9wl m w 
for weak decays in the Standard Model (SM). Consequently, the new flavor-changing 
couplings must be much smaller than the CKM mixing angles in order for this model to 
be phenomenologically viable. The most stringent bounds arise from the radiative decay 
B — > X s 7, which we discuss in the present work. (Contributions of light b squarks to 
kaon decays, K-K mixing, and D-D mixing are strongly suppressed.) The presence of 
such tight bounds implies stringent constraints on model building. 

If the light b squark is sufficiently light to be pair-produced in b decays, new uncon- 
ventional decay channels would be opened up, which could affect the phenomenology of 
B mesons and beauty baryons. Examples of potentially interesting consequences include 
modifications of beauty lifetime ratios, an enhancement of the semileptonic branching 
ratio of B mesons via production of charmless final states containing b squarks, an en- 
hancement of AT(B) and of the semileptonic CP asymmetry A^l, and wrong-sign kaon 
production via b — > sbb transitions allowed by the Majorana nature of the gluino. The 
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phenomenology of such effects will be discussed elsewhere. If the light b squark is too 
heavy to be pair-produced, it would still give rise to potentially large virtual effects in 
B decays. Their study is the main purpose of this Letter. 

2 The low-energy effective AB = 1 Hamiltonian 

We denote by di with i = 1, ..,6 the down-squark mass eigenstates, and by and q R 
with q = d,s,b the interaction eigenstates (the superpartners of the left-handed and 
right-handed down quarks). They are related by a unitary transformation q~L = Y } di 
and q~R = Vj di. We identify J3 with the light sbottom mass eigenstate and define a 
sbottom-sector mixing angle 9 through = cos# and r^ 3 = sin#. The fact that the 
light b squarks are not produced in Z decays implies 

sin^ ± yf si n^ ± ^(l-J). (1) 

The phenomenologically favored range for the Zd 3 d 3 coupling is | sin 6* | = 0.3-0.45 |J, 
meaning that the light sbottom is predominantly the superpartner of the right-handed 
bottom quark. In our numerical analysis we will assume a vanishing tree-level coupling 
to the Z and thus use sm9 = ±0.395. 

The flavor-changing couplings involving the light b and g fields can be parameterized 
by dimensionless quantities 

4 B = r S rg , (with ef = ef tan 6) (2) 

where A,B = L, R, and q = s or d for b — > s or b — > d transitions, respectively. In 
general the parameters e qh B are complex, which can lead to new CP-violating effects. 
These parameters are invariant under a phase redefinition of the light 6-squark state, 
and they transform in the same way as the products V* q Vib (with i = u,c, t) of CKM 
matrix elements under a phase redefinition of the down-type quark fields. It follows that 
ratios of the type e^/e^p and £f^/(V* q Vib) are invariant under phase redefinitions, and 
thus can carry an observable, CP-violating phase. 

Flavor-changing hadronic processes in the model with a light gluino and a very light 
b squark are most transparently described by means of an effective "weak" Hamiltonian. 
If we neglect effects that are suppressed by inverse powers of the heavy SUSY scale, the 
relevant energy scales are the electroweak scale, at which the usual SM flavor-changing 
operators are generated by integrating out the top quark and the W and Z bosons, and 
the scale m§, at which new flavor- changing operators are generated by integrating out the 
gluinos. We start by discussing the construction of the effective theory below the gluino 
scale, focusing on the new interactions proportional to e qb B induced by gluino exchange, as 
illustrated in Figure ffl. SUSY modifications of the renormalization-group (RG) evolution 
of the standard weak-interaction operators will be discussed later. The remaining light 
degrees of freedom in the low-energy theory are the quarks u, d, s, c, b, the photon and 
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b g s 



b s 




Figure 1: Examples of b — > s transitions induced by gluino exchange. The diagrams 
on the right show the corresponding contributions in the effective theory where the 
gluinos are integrated out. 



gluons, and the light b squark. Operators in the effective Hamiltonian can be organized in 
an expansion in inverse powers of the gluino mass. For m g « 15 GeV ^> rrib, it is a good 
approximation to keep the leading terms in this expansion, which have mass dimension 
five. These operators comprise the usual electromagnetic and chromomagnetic dipole 
operators, and new operators containing two scalar b fields. The effective Hamiltonian 
for b — > s transitions is (here and below, m g = rrig(m,g) denotes the running gluino mass 
at the gluino matching scale) 



TT 

a. 



susy _ ^a s {m~ g ) 



cff 



m- g 



\J*<% r <jjl) + ( L ^ R )} + O(VmJ) , 



(3) 



where 



0^ R = s L t a bb*t a b R , 



LR 



16tt 2 



s L cTuuF^b R , 



^ = s L bb*b R , 

LR 9s 
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16tt 2 
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(4) 



We work with the covariant derivative iD^ = id^ + eQdA 11 + g s A^t a , where Qa = — | is 
the electric charge of a down- type (s) quark. (To facilitate comparison with the literature, 
which usually adopts the opposite sign convention for the couplings, we have included a 
factor of —1 in the definition of the dipole operators O-j and Og.) In addition, there are 
dimension-five fermion-number violating interactions of the form s c (l ±75)66*6*, which 
mediate b — > sbb transitions. They are irrelevant to our discussion here. 

The Wilson coefficients at a scale fx ~ m g are obtained by matching the effective 
theory to the full theory. At leading order we find 



dirrig) = 2 , C 2 (m- g ) = , C 7 (m- g ) = Q d 
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4N 
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(5) 
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Figure 2: Two-loop diagrams relevant to the mixing of into 07 ;8 , and corre- 
sponding results, in units of o; s /47r, after the subtraction of sub divergences. Mirror- 
symmetric graphs with the gluon attached to the s-quark line give identical contri- 
butions. Results in the first column of the table have to be multiplied by the color 
and charge factors in the remaining columns. 



where N = 3 is the number of colors. In order to use the effective Hamiltonian for 
calculating B-decay amplitudes, we compute the values of the Wilson coefficients at a 
low scale /1 ~ m& by solving the RG equation (d/dln/i — 7 T ) C(fi) = 0. At leading order, 
the anomalous dimension matrix 7 receives contributions from the one-loop mixing of the 
operators (Oi, 2 ) and (0 7 , 8 ) among themselves, and from the two-loop mixing of Oi 2 
into 07,8- This is analogous to the case of the SM, in which one needs to consider the two- 
loop mixing of the current- current operators into the dipole operators at leading order 
fllfl . In our case, only the three two-loop diagrams shown in Figure |2| give a nonvanishing 
contribution. All other graphs vanish after their sub divergences are removed. The 
calculation of the UV divergences of these diagrams can be reduced to the evaluation of 
massive tadpole integrals |12|. The resulting anomalous dimension matrix in the operator 
basis (Oi, 2 , 7 , 8 ) reads 



/-6N+j^ -f + 2^ 2 " Qdii-^i) 
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47T 
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N — — 
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+ 0(a 2 s 



(6) 



The scale dependence of the Wilson coefficients is now readily obtained by solving 
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the RG equation. Setting N = 3, we find 



CM = y T 8 + jj v- 7/2 , cm = 1 (rr 8 - rr 7/2 ) , 

n / N ^ Z' 4 S 4 7/2 4 38 2/ , 1224 4/ o 
C 8 (p) = ( 1 ^ - 1 _ H£ 2/3^) (?) 

V39 60 ' 260 ' J K 1 

Here 77 = [a s (mg)/a s (/i)] 1//3o(5A1) , and 

x 11Ar 2 2iV 1 
f3o{n f ,n g ,n s ) = — n f - —n g - - n s (8) 

is the first coefficient of the generalized QCD (3 function in the presence of nj light Dirac 
fermions, n g light gluino octets, and n s light complex scalars. Numerical results for the 
coefficients Ci(fi) will be given in Table [l] below. The scale dependence of the Wilson 
coefficients below the scale rrig arises mainly from the mixing of 0\ with O2 and O7 with 
08- The mixing of 0% and O2 into the dipole operators turns out to be small numerically. 

The presence of light SUSY particles also affects the RG evolution of the SM con- 
tributions to the effective weak Hamiltonian below the electroweak scale. We will now 
discuss these effects for the operators of relevance to radiative B decays. 



3 The radiative decay B — ► X s ~f 

The inclusive radiative decay B —>■ is one of the most sensitive probes of physics 
beyond the SM. Indeed, we will see that this decay provides very stringent bounds on the 
flavor- changing couplings ef b H and ef^. The SM prediction for the B — > X s ^ decay rate 



is known at NLO [|T^, [14], 15, 16] and, within errors, agrees with the data. The change 



of this prediction due to the light SUSY particles present in our model is fourfold: 

1. The main effects are the genuine SUSY flavor- changing interactions due to quark- 
squark-gluino couplings. For /x < m§, these interactions are described by the 
effective Hamiltonian constructed in the previous section. 

2. Even in the absence of flavor- changing couplings in the SUSY sector (i.e., for 
vanishing e4, B ), the SM operator basis gets enlarged by dimension-six penguin 
operators with field content sbb*iD and sbgg. These new operators have small 
Wilson coefficients and yield a negligible contribution to the B — > X s 7 decay rate. 
In our analysis, we will neglect these as well as all four-quark penguin operators. 

3. The presence of a light gluino octet (and, to lesser extent, of a light b scalar) mod- 
ifies the running of the strong coupling constant. We use the two- loop expression 
for a s (fi), modified to account for the effects of the light SUSY particles. 
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4. There is a SUSY contribution to the 6-quark wave-function renormalization, which 
adds to the anomalous dimensions of the SM operators. Also, the masses of the b 
quark and the gluino mix under renormalization, thus altering the scale dependence 
of rrib(fi). This is a novel effect due to the decoupling of the heavy bn squark at 
the SUSY scale. 



The last two effects change the evolution of the Wilson coefficients of the SM operators 
between the electroweak scale and the scale /i ~ m§, where the gluino degrees of freedom 
are integrated out. (Beyond leading order, the anomalous dimensions of the SM operators 
are also changed due to internal loops involving SUSY particles.) 

The effective weak Hamiltonian governing B — > X s 7 decays in the SM is 

AC 

H^ k = v*v tb CM QM ■ (9) 

The operators relevant to our calculation are 

Qx = AlA Al^L , Q2 = s Ll ^ L c L ^b Ll Q 7l = m b O LR , Q 8g = m b O LR . (10) 

To an excellent approximation, the contributions of other operators can be neglected. 
To obtain the values of the corresponding Wilson coefficients Ci(fi b ) in our model, we 
first evolve them from the electroweak scale \i = my/ down to a scale fi g ~ rrig, and in a 
second step from \x g to a scale fib ~ ra b - 

Above the gluino scale, there are SUSY contributions to the wave-function renormal- 
ization constants of left- and right-handed 6-quark fields from gluino-squark loops. At 
one-loop order, we obtain the gauge-independent results 

5Z 2 (b L ) = sin 2 £, 8Z 2 (b R ) = -°f± cos 2 6. (11) 

Next, by calculating the self-energies of b quarks and gluinos we find that their masses 
mix under renormalization. The corresponding anomalous dimension matrix in the basis 
(mb,m g ), defined such that drrii/dhifj, = —('J m )ij' m ji rea ds 



a. 
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Note that the off-diagonal entries are sensitive to the sign of the mixing angle 9. The 
RG evolution of the operators Qi in @ is complicated by the effects of mass mixing. To 
compute the resulting modifications of the Wilson coefficients we work in the extended 
operator basis (Qi,Q2,mbOj R ,m b Og R ,m g 7 R ,rrigOQ R ). Using (|TTJ), (0)> an d the well- 
known anomalous dimensions of the SM operators Jll|], we obtain for the anomalous 
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dimension matrix 7g = f^^n + O(o%), where (setting N = 3 and Qd = —\] 

/§sin 2 fl-2 6 3 \ 
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sin 29 
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6 §sin 2 £-2 ^ f 7 

f - I sin 2 9 I sin 29 



(o) 

7(i -f f-fsin 2 fl | sin 29 

sin2# f-fsin 2 # 

sin2# -f ±fi-fsin 2 #/ 

(13) 

The solution of the RG equation in this basis yields coefficients (c\, 02,0^,04,0^,0%) at 
a scale between and m g . Their initial values at the electroweak scale are given 
by (0, 1, C 7l (m w ), Cg g (mw), 0, 0). The relevant /3-function coefficient in this range is 
/5 (5, 1, 1). From these solutions, we obtain the SM Wilson coefficients at the scale 
fig ~ rrig by means of the relations Ci^(fig) = Ci j2 (/ig) and 

CV>,) = c 3 (M§) + ^4 C5 ^) , = c^g) + ^4 ■ (14) 

m b {Hg) m b {Hg) 

The sign of the coefficients 05^ depends on the sign of the mixing angle 9. At lead- 
ing order, the running 6-quark mass at the gluino scale is obtained from m^/i^) = 
m.K)^^)/^^)] 4 /*^ ' 1 ). 

Once we have determined the SM contributions to the Wilson coefficients at the 
scale fig, their evolution down to lower scales is governed by the well-known evolution 
equations of the SM. The corresponding 4x4 anomalous dimension matrix coincides 
with the upper left 4x4 corner of the extended matrix in (|13"D evaluated at 9 = 0. 
The resulting formulae are more complicated than in the SM, because in our case the 
coefficient Ci(fi g ) does not vanish at the matching scale (whereas Ci(m w ) = for the 
standard evolution). 

Table [I] shows the results for the Wilson coefficients at different values of the renor- 
malization scale. (The coefficient C\ does not enter the B — > X s 7 branching ratio and 
is omitted here.) The values of C 7l and C$ g depend on the sign of the mixing angle 9, 
although this effect is numerically small. For comparison, the values obtained at /1 = m;, 
in the SM (using a s (m z ) = 0.118) are C 2 ^ 1-12, CV 7 ~ -0.32 and C 8g ~ -0.15. 
They are very close to the values found in the presence of the light SUSY particles. In 
addition, there are the extra contributions proportional to the new coefficients Cj. The 
second column in the table shows the running 6-quark mass at the various scales. Note 
that the value of above the gluino scale is very sensitive to the sign of 9. The result 
rrib{mw) = 2.75 GeV corresponding to positive 9 appears to be favored by the DELPHI 
measurement m b {mw) = (2.67 ± 0.50) GeV obtained from three-jet production of heavy 
quarks at LEP ||17|| . However, here a similar comment as in our discussion of the run- 
ning of a s applies, namely that the DELPHI analysis implicitly assumes that m b runs 
as predicted in the SM. 
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Table 1: Results for the Wilson coefficients and the running 6-quark mass for 
different values of //. Input parameters are m b {m b ) = 4.2 GeV, m t (mw) = 174 GeV, 
m g{ m g) = 15 GeV, and a s (m b ) = 0.21. In the upper portion of the table the gluino 
is integrated out at /i — m g , in the lower portion at /i — 2m g . If two signs are 
shown, the upper (lower) one refers to positive (negative) mixing angle 9. 



Scale 


m b (n) [GeV] 


Ci 


c 2 


c 7 


c 8 


c 2 


C77 


C'sg 


mw 


3.17 T 0.42 










1 


-0.195 


-0.097 


m- g 


3.59 


2 





-0.222 


-0.833 


1.040 


-0.255 ±0.028 


-0.124 ±0.014 


m b 


4.20 


2.691 


0.066 


-0.264 


-0.804 


1.104 


-0.313 ±0.023 


-0.143 ± 0.011 


m w 


3.13 T 0.25 










1 


-0.195 


-0.097 


m g 


3.59 


2.274 


0.025 


-0.241 


-0.821 


1.042 


-0.255 ±0.015 


-0.124 ±0.008 


m b 


4.20 


3.064 


0.104 


-0.280 


-0.791 


1.106 


-0.313 ±0.013 


-0.143 ± 0.006 



We are now ready to present our results for the B — > X s ^ decay rate, including both 
the SM and the new SUSY flavor-changing contributions. It is convenient to define new 
coefficients 

r i>R(,,\-r (,.\ ^ a s{mg) e L sb R C 7 (/x) 

c 7 w-g,m- Gpm _ 

67 v^MiT)' (15) 

and analogous coefficients C% R and C$ L . These expressions exhibit the general fea- 
tures of our model as described earlier. The SUSY contributions are enhanced relative 
to the SM contributions by a large factor y/ ^ na ^ m ^ ~ xq 3 , meaning that the ratio of 

flavor-changing couplings, e^, B /(U t * V tb ), must be highly suppressed so as not to spoil the 
successful SM prediction for the branching ratio. The resulting leading-order expression 
for the B — * X s ^ decay rate is 



G 2 F aM*m 2 b {m b ) 



T(B - X sl ) = t ^-p 21 I^I^IC^^)! 2 + |C 7 ^M| 2 j , (16) 

where \i b ~ m b is the renormalization scale. M b is a low-scale subtracted quark mass, 
which naturally enters the theoretical description of inclusive B decays once the pole 
mass is eliminated so as to avoid bad higher-order perturbative behavior. 

It is well known that NLO corrections have a significant impact on the B X s 7 
decay rate in the SM, which is largely due to NLO corrections to the matrix elements 
of the operators Qi in the effective weak Hamiltonian. (NLO corrections to the Wilson 
coefficient C 77 have a much smaller effect.) In order to capture the bulk of these correc- 
tions, we include the 0(a s ) contributions to the matrix elements but neglect SUSY NLO 
corrections to the coefficients C? R and C RL . We also neglect two-loop contributions to 
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Table 2: Results for the coefficients Bo and Ai j2 for the SM (first row), and for 
the SUSY scenarios with positive (middle portion) and negative (lower portion) 
mixing angle 6, for E 1 > 2 GeV. The quoted errors refer to the variations of the 
theoretical parameters within the ranges specified in the text. The renormalization 
scale is varied between 2.5 and 7.5 GeV. Other input parameters are V t * s V t b = —0.04, 
Oi s (mb) = 0.21 for the SUSY scenario, and a s (mb) = 0.225 for the SM. 





Default 


A((ml s ) 3 m 2 b ) 


A(m c /m b ) 


Afi b 


c>SM 
n 


3.44 


±0.21 


TO. 11 


+0.09 
-0.16 


B 


2.93 


±0.18 


T0.09 


+0.06 
-0.13 


10- 4 Ai 


3.60 





±0.05 


+0.12 
-0.02 


10- 8 A 2 


3.12 





±0.10 


+0.28 
-0.12 


Bo 


3.68 


±0.23 


TO. 11 


+0.04 
-0.11 




3.19 





±0.04 


+0.06 
-0.01 


1Q~ 8 A 2 


2.48 





±0.07 


+0.18 
-0.08 



the matrix elements involving 6-squark loops. This is justified because of the relatively 
large mass of the b squark, and because our two-loop anomalous dimension calculation 
has shown that there is very little mixing of the squark operators into the dipole opera- 
tors. At NLO our results become sensitive to the precise definition of the mass parameter 
Mb, which we identify with the so-called Upsilon mass ||18|| , for which we take the value 
ml = 4.72 ± 0.06 GeV (Up to a small nonperturbative contribution, m\ s is one 
half of the mass of the T(IS') resonance.) We also introduce a cutoff E™ m = — S) ml 
on the photon energy in the 5-meson rest frame, which is required in the experimental 
analysis of radiative B decays. We then obtain 



Br(5^ X sl )\ 



G F a{m b 



1S\3 



ml(m b ) 



32tH 



v; s v tb \ 



(17) 



where K^q(S) is obtained from the formulae in [ff5| by obvious modifications to include 
the effects of the new SUSY contributions, and by a change in some of the NLO terms due 
to the introduction of the Upsilon mass in flT6| ) pOfl . The dependence of the branching 
ratio on the SUSY flavor- changing couplings can be made explicit by writing 



Br(B - X sl ) 



Ey>E" 



10- A B (5) 1 T A 1 (5) Re(e^) T A 2 {8) 



le^l 2 + 



ML\2 



(18) 

In Table 0, we give results for the coefficients Bq and A\ } 2 including the dominant theo- 



retical uncertainties. Following |T6[ , we use a running charm-quark mass in the penguin- 
loop diagrams rather than the pole mass. This is justified, because the photon-energy 
cut imposed in the experimental analysis prevents the intermediate charm-quark prop- 
agators from being near their mass shell. Specifically, we work with the mass ratio 
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Figure 3: Allowed regions (at 95% c.l.) for the SUSY flavor- changing parameters 
obtained from the CLEO measurements of the B — > V s7 (left) and B X sg (right) 
branching ratios, using central values for all theory input parameters. The shaded 
regions correspond to the SUSY model with positive mixing angle 9, the dashed 
lines refer to negative 9. 



m c (fj,) / rrib(fi) , where the running masses are obtained from m c (m c ) = (1.25 ± 0.10) GeV 
and m b (m b ) = (4.20 ± 0.05) GeV. 

In the left-hand plot in Figure |3|, we confront our theoretical result for the B — > X s j 
branching ratio with the CLEO measurement By(B — ► X s j) = (3.06 ±0.41 ±0.26) ■ 10~ 4 
obtained for E 7 > 2 GeV |21] . (This result actually corresponds to the sum of B — > X s j 
and B — > X^7 decays. However, the suppression of the exclusive B — > p 7 decay with 
respect to the B — > K*^ mode implies that the dominant contribution to the inclusive 
decay must come from b — > s 7 transitions. In the context of our model, it follows that 
the couplings ej[ B must obey even tighter constraints than the ef B .) It follows that the 
maximum allowed values of the parameters ef B are 10~ 4 , as is already obvious from the 
magnitude of the coefficients Ai t % in Table [| However, values larger than 5 ■ 10~ 5 would 
require a fine-tuning of the phase of e^ b R and are thus somewhat unnatural. Note that 
the ratio shown on the vertical axis in the plot is bound to lie between 1 and —1, and in 
the limit e B b = corresponds to cos$lr, where denotes the CP- violating phase of 

e sb ■ 

The right-hand plot in the figure shows a similar constraint arising from the inclusive 
charmless decay B — > X sg . At leading order, the decay rate for this process is obtained 
from (0) by the replacements a — > ^ot s {ii b ) and Cj B — > C% B . The allowed region 
corresponds to the CLEO upper bound of 8.2% (at 95% c.l.) for the B — > X sg branching 
ratio [P^j . There are, however, potentially large theoretical uncertainties in this result, 
because we neglect NLO corrections to the branching ratio. We therefore refrain from 
combining the two plots to reduce the allowed parameter space. 
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In the SM, the direct CP asymmetry in the inclusive decay B — > X s 7 is very small, 
below 1% in magnitude p3| . In the SUSY scenario, on the other hand, the phase of the 



coupling could lead to a large asymmetry. In the approximation where one neglects 
the SUSY contributions to the CP-averaged decay rate in (|16D , which is justified in view 



of the good agreement of the SM prediction with the data, the formulae in [23] yield the 
prediction Aqp w —50% x 10 4 Im(ef fe ' R ), where we have neglected the small contribution 
from the charm-quark loops and the yet smaller contribution from 6-squark loops. (We 
use the standard phase convention where X t = V t * V t b is real. In general, the CP asym- 
metry depends on Im(e^ 6 /At).) It follows that even within the very restrictive bounds 
shown in Figure [] there can be a potentially large contribution to the CP asymmetry, 
which would provide a striking manifestation of physics beyond the SM. In fact, the 
CLEO bounds -27% < A CP < +10% (at 90% c.l.) g| imply that 

-2 • 10" 5 < Im(eff) < 5 • 10" 5 , (19) 

placing another tight constraint on the flavor-changing coupling e^. 



4 Conclusions 

New supersymmetric contributions to 6-quark production at hadron colliders can ac- 
count for the long-standing discrepancy between the measured cross sections and QCD 
predictions if there is a light b squark with mass in the range 2-5.5 GeV, accompanied 
by a somewhat heavier gluino 0. In this Letter, we have explored the phenomenology 
of rare B decays in such a scenario and have found tight bounds on the flavor-changing 
parameters controlling supersymmetric contributions to b — > s and b —>■ d FCNC tran- 
sitions. The most restrictive constraints arise from virtual effects of light b squarks in 
B — > X s 7 decays. We have analysed this process by constructing a low-energy effective 
Hamiltonian, in which the gluinos are integrated out, while the b squarks remain as 
dynamical degrees of freedom. We find that the flavor- changing couplings ef b L and ef b R 
must be of order few times 10~ 5 or less. (Even tighter constraints hold for the analogous 
b — > d couplings.) This implies that certain off-diagonal entries of the down-squark mass 
matrix must be suppressed by a similar factor compared to the generic squark-mass 
squared. 

Even with such tight constraints on the couplings, this model allows for interesting 
and novel New Physics effects in weak decays of B mesons and beauty baryons. As 
an example, we have discussed the direct CP asymmetry in B —>■ X s 7 decays, which 
could be enhanced with respect to its Standard Model value by an order of magnitude. 
Other possible effects include an enhanced B — > X sg decay rate. We have not considered 
here the possibility of 6-squark pair-production, which would be kinematically allowed 
for very light squark masses. The new decay modes b — *> sbb* and b — *> sbb would affect 
the decay widths of B mesons and A& baryons differently, and hence might explain the 
anomaly of the low A& lifetime. We will report on this interesting possibility elsewhere. 
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